FREQUENCY SHIFTS AND LINEWIDTH CHANGES OF 
INFRARED-ACTIVE PHONONS IN DOUBLE-LAYERED 
HIGH-TEMPERATURE SUPERCONDUCTORS 
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We calculate frequency shifts and changes in linewidths of infrared-active phonons 
within a shell model for the bare phononic system coupled to an electronic double- 
layer structure with inter-layer charge transfer. The theoretical concept is applied 
to YBa 2 Cu307 yielding a good description of experimental results in the normal 
state as well as at the transition to the superconducting state. 
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1 Introduction 



Raman- [[IJ |2| and infrared-experiments || || show 
a strong influence of the superconducting transition on 
phonon frequencies and linewidths in YBa 2 Cu 3 07 and 
similar compounds. The experimental results on Raman- 
active phonons are explained by the theory of Zeyher 
and Zwicknagl || |Kj . Complementary to their work the 
even greater variety of effects found on infrared-active 
phonons, e. g. narrowing of phonon-lines [[jj and soft- 
ening of high-frequency phonons || - || in connection 
with the transition into the superconducting state, has 
been explained qualitatively [11], [13] within a double- 
layer model for high-temperature superconductors. 
In this paper we link the approaches of refs. [12| , where 
a microscopically based model for a Coulomb-type elec- 
tron-phonon interaction in layer systems was developed, 
with the calculation of the electronic susceptibility for 



the double-layer model from ||13|| . This results in a mi- 
croscopically founded description of phonon line width 
changes and frequency shifts which fits to the experi- 
mental data in the normal state as well as in the super- 
conducting state. 



2 Infrared absorption in coupled 
phonon-layer systems 

We are going to investigate the infrared absorption of an 
essentially ionic system coupled to conduction electrons 
which are confined to two-dimensional planes within a 
model presented in [12, |13| : Bare phonons are described 



by a shell-model of lattice dynamics, bare electrons are 
allowed to move almost freely within the conducting 
Cu02-planes; we assume hopping in c-direction between 
different planes within one unit cell but neglect trans- 
port between different double-layers for simplicity. We 
consider electron-phonon coupling due to the Coulomb- 
interaction between vibrating ions and conduction elec- 
trons. 

In order to determine the response function for IR ab- 
sorption, the dielectric function, of a coupled electron- 
phonon system it is not sufficient to calculate a renor- 
malized phononic Green's function and, independently, 
an — also renormalized — electronic susceptibility. E- 
lectron-phonon coupling yields "off-diagonal" terms in 



the dielectric function which may change lineshapes of 
absorption lines. 

This can be described by a generalized susceptibility of 
the coupled electron-phonon system depending on fre- 
quency u and wavevector q. Its inverse is given by 



X£(q,") = ( G * J-x ] ■ (1) 



Here, G p h is the bare phonon Green's function which is 
determined by 

l(u + i5) 2 - D ] M 1/2 G ph M 1/2 = 1 (2) 

with the dynamical Matrix Do(q) and diagonal matrices 
M containing the vibrating masses. For reasons of sim- 
plicity we do not make any formal distinction between 
cores and shells in the lattice dynamical shell-model used 
for the description of the bare phonon system. Instead, 
we use shell masses about 10 6 times smaller than core 
masses. Thus, for a system of N ions Do is a (QN x QN) 
matrix whose 3N lowest eigenvalues describe the disper- 
sion of the bare vibrational system. 
x(q, uj) is the susceptibility of the bare electronic sys- 
tem. In general this is a (n x n) matrix where n is the 
number of conducting layers per unit cell. In the case 
of a double-layer system and q parallel to the planes 
this matrix can be transformed into diagonal form with 
elements describing intra- and interband excitations, re- 



spectively |13 



Finally, 7(q) in (H) is a (n x 6N) matrix containing the 
electron-phonon matrix elements. We calculate these 
matrix elements according to our ionic approach de- 
scribed in detail in [|l^] (with X MK from [[12] and ionic 
charge numbers Z K we have -y^ K = JB^X^ for fj, = 1,2, 
k = 1,...,6N). 

In order to calculate the infrared absorption we still have 
to multiply the matrix elements of Xgen by the coupling 
constants to the external electric field. For the case of 
infrared active phonons with eigenvectors perpendicu- 
lar to the a-6-plane and charge fluctuations between the 
double layers the dielectric function reads 

Ae(u) = 4vrA + x g en(q -> 0, u)A (3) 

with 



A K = eZ K fur k = 3v; v = 1, . . . 2N (4) 



Aqn + i = ±0.5ed 
Aqn+2 = =F0.5erf 



(5) 
(6) 



All other components of the external field coupling vec- 
tors A are zero. The electrical field couples to the polar- 
ization in the crystal. For lattice vibrations this polar- 
ization is given by the product of ionic charge (or shell 
charge) and the corresponding displacement. The latter 
is contained in x gen - For the electronic excitations the 
polarization is determined by the distance d between the 
two planes connected by the electronic hopping process 
and the charge transfered between these layers which 
depends on the electronic susceptibility. The signs and 
distances to be used in (|5|) and (|6|) depend on the way the 
inter-layer hopping process is taking place physically. 



3 Infrared absorption of YBa2Cu30^ 

In the following we apply the theoretical concepts de- 
scribed above to YBa 2 Cu 3 7 using a lattice-dynamical 



shell model with parameters from |14|| . As we are mainly 
interested in general trends we refrain from a "fine-tun- 
ing" of the original shell-model parameters. 
The electronic double-layer system is modelled by two 
bands crossing the Fermi surface Jl3| . For simplicity the 



bands are assumed to have zylindrical symmetry but an 
energy distance varying from zero to about 1000 cm -1 
in the vicinity of the Fermi surface. The mean density of 
states is 0.85 eV _1 for each band. In the superconduct- 
ing state we use a BCS-type approximation with equal 
gaps 2A = 320 cm" 1 in both bands. 
The bare electronic susceptibility is calculated numer- 
ically as described in |L3] taking into account vertex 
corrections and screening effects in a random-phase ap- 
proximation. We assume that the electronic hopping 
process takes place via the copper-oxygen chains, while 
the yttrium layer is considered to be completely insulat- 
ing. Thus, the distance d to be used in (H) and © is 
about 8.3 rA. 

The electron-phonon coupling at large wavelengths is 
completely determined by the Coulomb interaction be- 
tween ionic and electronic charges without additional 
parameters as described in [l^ . 



Fig. la shows the result of the calculation for the nor- 
mal state. The dashed line represents the shell model 



calculation (with an artificial broadening of the phonon 
lines), the solid line represents the results for the cou- 
pled system including the electron-phonon interaction. 
Three prominent effects are to be observed: softening 
of phonon frequencies, broadening of absorption lines 
and strong changes of oscillator strengths for different 
modes. 

Most pronounced are these effects for the mode at 365 
cm -1 in the uncoupled system, where the yttrium-layer 
moves against the conducting CuCVplanes flH| . The 
frequency of this mode is lowered by about 85 cm -1 due 
to the electron-phonon coupling, the linewidth increases 
as well as the oscillator strength. Actually, such a fre- 
quency shift for this mode is observed experimentally 
going from YBa2Cus06 to YBa 2 Cus07 Similar ef- 

fects are found also for the modes at 154 cm -1 and 508 
cm -1 (in the bare model) which both become softer by 
13 cm -1 . 

These softenings result from the fact that ionic displace- 
ments perpendicular to the conducting double layer in- 
duce a charge transfer between the planes of this double 
layer. This charge transfer is accompanied by a polar- 
ization field |12| which reduces the restoring forces onto 
the ion, the frequency of the vibration becomes lower. 
The line broadening is caused by the damping of lattice 
vibrations due to interband particle-hole excitations in 
the electronic system. 

Somewhat less evident is the situation concerning the 
changes in the oscillator strength. Here, one has to 
distinguish two cases: firstly, vibrations of ions which 
are located between the two planes connected by the 
electronic hopping process, in our case this would be e. 
g. the barium ions; here the direction of the electronic 
polarization due to the hopping process is antiparallel 
to the polarization resulting from ionic displacements, 
therefore the oscillator strength of the mode (around 150 
cm -1 ) decreases. If on the other hand e. g. yttrium ions 
move perpendicular to the planes, the induced electronic 
polarization is in the same direction as the phononic 
polarization, the oscillator strength therefore increases. 
Experiments also show large oscillator strengths for bar- 
ium modes || in contradiction to our results. This may 
be due to additional charge transfer processes, e. g. in 
the CuO-chains, which are not contained in our model. 
Our numerical results for the oscillator strength depend 
to a large extent on the details of the electronic band 



structure, while the results for the frequency shifts are 
more independent of such details. Our calculations for 
the dielectric function also yield a direct contribution 
from collective electronic charge fluctuations between 
the layers. It appears at much higher frequencies and 
therefore is not shown in Fig.l. 

What happens at the transition to the superconducting 
state! The infrared absorption spectrum in the super- 
conducting state is shown in Fig. lb. In Figs. 2a and b 
we show the relative changes in frequency and linewidth 
due to this transition calculated within the framework of 
our model. The small symbols refer to different experi- 
mental results, the circles are the results of our numeri- 
cal calculations and fit well the experimental frequency 
shifts. All phonons — even those at high frequencies 
- become softer in the superconducting state, phonon 
lines below 500 cm -1 become narrower due to the loss of 
decay channels resulting from the gap formation in the 
superconducting state. Above 500 cm -1 damping in- 
creases and phonon lines may become broader in the su- 
perconducting state. The strongest effects are observed 
for the modes arround 300 cm -1 . 

We have also calculated phonon spectra at finite wave 
vectors parallel to the planes. Here we find that the 
relative frequency shifts and line width changes at the 
transition to the superconducting state become much 
smaller. This is confirmed by recent neutron scattering 
experiments [TH . 



4 Summary 

In this work we have presented calculations for the in- 
frared absorption of YBa 2 Cu 3 07 using a microscopically 
based double-layer model for high-temperature super- 
conductors. 

The approximations used to describe the conduction 
electron system are rather crude in view of strong-corre- 
lation effects. Nevertheless, they allow for a treatment 
of realistic systems and we regard the results to be qual- 
itatively correct in the long wavelength limit. 
The electron-phonon coupling was described without any 
additional parameters by the Coulomb interaction be- 
tween the ions inside and outside the Cu0 2 -planes and 
the mobile electrons. Again, we believe that this type of 
interaction will dominate the interband electron-phonon 



coupling in the limit of small wavevectors. 
Our model calculations explain experimental results on 
frequency shifts and line width changes of infrared-active 
c-axis phonons. In particular the softening of modes at 
the transition to the superconducting state is well de- 
scribed by our model. Concerning the absolute postion 
of the absorption peaks the agreement with experimen- 
tal results could be improved if we would adjust the 
parameters of the original shell-model. We did not do 
this, because we think it more important to first improve 
the approximations for the electronic system. In some 
experiments |7[] it seems, that the softening of phonon 
frequencies sets in already at temperatures above the su- 
perconducting transition temperature. At the moment 
it is not clear wether this effect can be explained by su- 
perconducting fluctuations or by some other mechanism 
like opening of a spin gap. 
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Figure Captions 



Fig. 1. Calculated infrared absorption for YBa2Cus07 
without electron-phonon interaction (dashed line) and 
with electron-phonon interaction (solid line); a) results 
in the normal state (T=100K) and b) in the supercon- 
ducting state T=1K. 

Fig. 2. Frequency shifts (a) and line width changes (b) 
of infrared-active c-axis phonons in i?Ba2Cus07 (R = 
rare earth). The small symbols are experimental results 
from refs. [[| - || . The circles mark the frequency shifts 
and line width changes calculated in this work. 



